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Abstract: Very-high-energy (VHE, E > 100 GeV) 7 radiation has already been detected from several supernova remnants 
(SNRs). These objects, which are well-studied in radio, optical and X-ray wavelengths, constitute one of the most 
intriguing source classes in VHE astronomy. H.E.S.S., an array of four imaging atmospheric Cherenkov telescopes in 
Namibia, has recorded an extensive dataset of VHE 7-ray observations covering the central region of the Milky Way, 
both from pointed observations as well as from the Galactic Plane Survey conducted in the inner region of the Galaxy. 
From radio observations, several hundred SNRs are known in the Milky Way, but until now only few of them have been 
identified as VHE 7-ray emitters. Using the H.E.S.S. dataset and a large ensemble of radio SNRs localized in the inner 
region of the Galaxy, the standard framework that links the origin of cosmic rays to the 7-ray visibility of SNRs can now 
be tested. Here we present the ensemble of investigated SNRs and discuss constraints on the parameter space used within 
a theoretical model of hadronic VHE 7-ray production. 
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1 Introduction 

SNRs are widely believed to be the main sources of Galac- 
tic Cosmic Rays (CRs) with energies of up to 10^^ eV. The 
amount of energy released during a supernova (SN) explo- 
sion, the rate of Galactic SNe and the energy density of 
CRs together with diffusive shock acceleration as a pos- 
sible acceleration mechanism make SNRs the prime can- 
didates for the acceleration of CRs. In radio, optical and 
X-ray wavelengths deep observations on SNRs have been 
performed and several SNRs have been discovered in the 
relatively young field of VHE 7-ray astronomy . With 
H.E.S.S. the shell-type morphology of individual SNRs 
such as RX J1713.7-3946 1 1 1, SN 1006 [2|, Vela Junior [3| 
and HESS J1731-347 p-l has been resolved. These discov- 
eries support the hypothesis of particle acceleration in the 
shell of SNRs. Using the large dataset of H.E.S.S. accumu- 
lated during the Galactic plane survey (GPS) [5 1 in recent 
years we study SNRs in the inner Galactic region, ranging 
from 275° to 60° in Galactic longitude and from -3° to 3° 
in Galactic latitude. 

Here, the VHE 7-ray signal of known radio SNRs listed 
in Green's radio SNR catalog f6l is investigated. For each 
SNR within the GPS region, an upper limit on the VHE 
7-ray flux has been calculated and can be compared to the 
model predictions. In this work we test the commonly used 
model of the hadronic 7-ray emission expected from SNRs 



presented by Drury et al. ff\. They calculate the integral 
7-ray flux from hadronic interactions based on intrinsic SN 
parameters. VHE 7 rays result from the decay of neutral pi- 
ons generated in proton-proton interactions in the expand- 
ing shell of the SNR. For the first time, the high sensitivity 
of H.E.S.S. and the large dataset of the GPS allow us to 
derive upper limits that are low enough to test model pa- 
rameters. Using the upper limit on the VHE 7-ray flux 
represents a conservative test since an additional leptonic 
component might be present in the flux. 

2 Probing the hadronic model 

Assuming that CRs are accelerated and VHE 7 rays pro- 
duced in SNRs, the flux of VHE 7 rays depends on ex- 
plosion and acceleration parameters, the properties of the 
ambient medium, and the distance of the SNR. Drury et 
al. |7| calculate the expected integral flux of VHE 7 rays 
from SNRs generated from hadronic CR interactions above 
a given energy threshold E: 
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where Q is the efficiency of the particle acceleration, /sgN 
the total energy released during the supernova explosion, d 
the distance to the SNR and n the density of the interstellar 
medium surrounding the SNR. Typical values are = 0.1, 
-Bsn = 10^^ erg and n = 1 cm^'^ |7|. For this model 
the authors assume a power law with a spectral index of 

r = 2.1. 

3 SNR sample 

203 SNRs listed in Green's radio SNR catalog are located 
within the H.E.S.S. GPS region. For these SNRs an up- 
per Umit on the VHE 7-ray flux can be derived based on 
H.E.S.S. observations. A subsample used for the investi- 
gation of the hadronic model is chosen following theoret- 
ical considerations. The essential criterion for VHE 7-ray 
production is the evolutionary stage of the respective SNR. 
Except for the very young SNR Gl.9+0.3 with an age of 
~1 10 years [81, all SNRs with known age can be expected 
to have reached the evolutionary stage of the Sedov phase 
where efficient particle acceleration takes place. An upper 
age bound for the VHE 7-ray visibility of SNRs has been 
discussed in several publications. Drury et al. |7| give an 
age of ^10"* years where temperatures drop below ^10^ K 
and energy losses due to thermal cooling and particle es- 
cape from the shell due to magnetic instabilities become 
more important. Ptuskin & Zirakashvili |9| study nonlin- 
ear effects in SNRs leading to a decrease of the maximum 
energy of accelerated particles. A generic prediction of 
the models they consider is a reduced VHE 7-ray visibil- 
ity of SNRs at an age larger than a few thousand years. In 
the present work we therefore restrict the investigation of 
the hadronic model to SNRs with a known age of up to 
5000 years. 

Predicting the VHE 7-ray flux of SNRs requires a reliable 
estimate of their distance. These distance estimates are 
usually obtained through the comparison of HI absorption 
spectra with respective HI and molecular (^^CO, ^'^CO) 
emission spectra along the line of sight (see [ lO] and ref- 
erences therein). The uncertainties of the Galactic rotation 
curve model used to estimate these kinematic distances, to- 
gether with those on the distance to the Galactic Center and 
the circular velocity of the Sun lead to errors of the order of 
10-30%. For relatively nearby (a few kpc away) SNRs with 
detected radio, optical or X-ray emission, it is sometimes 
possible to measure the expansion rate of the shell (or oth- 
erwise, that of the associated pulsar wind nebula, see 1 1 1 1 
for G2 1.5-0. 9) which constrains the dynamical age (subject 
to the evolutionary stage of the SNR) and, therefore, pro- 
vides information about the distance to the source (see [12] 
for Vela Junior). 

The 20 SNRs considered within this study are listed in Ta- 
ble[T]together with their mean radius R, distance d, age, the 
VHE 7-ray flux upper limit (<i>'^^, seed]) and the derived 
upper limit on the model parameters {9Ea,^vF^ , sec. [5]). 
The mean radius R is derived from the major and minor 
diameter given in Green's catalog as the geometric mean 



R = Y2-^I?niaj • £'min- The SNR positions have also 
been taken from the catalog. For the analysis the same 
spectral index of F = 2.1 has been adopted from the 
model of Drury et al. for all SNRs, allowing us to per- 
form a standardized analysis of the entire SNR sample. For 
some sources the spectral index is known to be other than 
2.1 which leads to differing analysis results compared to 
a detailed spectral analysis (see Aharonian et al. I.13J for 
G0.9H-0.1, HESS J1747-281). 

4 Upper limit on the VHE 7-ray flux 

For the calculation of a flux upper limit, a reference flux 
with a predefined spectral shape has to be assumed, 

c^,,i{E) = 00 {E/Eo)-^ (2) 

The normalization can be chosen arbitrarily as it will cancel 
out in the end. We use E^ = 1 TeV in the following. For 
any given run with live-time T, the number of expected 
signal events from a source can then be calculated as 

/>oo 

N,un = T / ^,,t{E)A,s{E,q)dE (3) 

Here, E is energy, A^s is the effective area, q symbolizes 
the observation parameters (e.g. the zenith angle) and E'min 
is the threshold energy appropriate for the run. The indi- 
vidual expected counts can then be summed over all runs 
to yield the total number N of expected counts. 

An upper limit on the integral flux above energy Ei is then 
calculated as 

UL roo 

$UL = — / cl>,,f{E)dE (4) 

J El 

where n^^ is the upper limit on the excess counts above 
threshold that can be obtained from the measured number 
of photon candidates above threshold within the signal re- 
gion, A^oN, the corresponding number of counts from a 
suitable background control region, -/Vqff, and the back- 
ground normalization factor a llT4ll . by using the profile 
likelihood method of Rolke et al. IfTSl . 

The signal region is taken to be a circle of radius 9, where 
e = 0.1° for SN Rs with radius R < 0.05° and 6 = 
VR^ + 3 • 6ipsF^ for SNRs with radius R > 0.05°. This 
means that small SNRs are treated as pointlike sources and 
6 is chosen such that the significance gain a/^/i is opti- 
mized for the point spread function (PSF) of the instru- 
ment. The PSF is characterized by the scale 6'psf — 0.1° 
here. For SNRs that have to be considered as extended on 
the scale of the PSF, the integration radius is chosen such 
that the emission of the SNR can be assumed to be fully 
contained and spillout from PSF effects is negligible. The 
effective area curves were calculated separately for the two 
cases. The flux upper limits are calculated at the 95% confi- 
dence level and given in units of the flux of the Crab nebula 
of $crab(i^ > ITeV) 2.2 • lO-^cm-^s"!. 
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5 Discussion 

Supernova remnants are the most plausible sites of cosmic- 
ray acceleration. As a consequence, they are predicted 
to be emitters of VHE 7 radiation by tt" decay during 
the early stage of their evolution. The good sensitivity of 
H.E.S.S. and the large dataset accumulated in its Galactic 
Plane Survey now allow us to compare quantitative esti- 
mates for this emission to observations, and to study the 
SNR population in the Inner Galaxy as a whole. From 
Green's catalog of radio SNRs, candidates were selected 
whose estimated age makes them predicted sources of VHE 
7 radiation and the existance of a distance measurement al- 
lows for the expected flux to be calculated. 

Fig. [T] shows the distribution of the upper limits on the 
product p = OEsj^sn/lO^^ ergcm^'^ as calculated from the 
H.E.S.S. flux upper limit in the framework of the model of 
Drury et al. The vertical hne marks a value of p — 0.1 
which would be expected for a supernova explosion en- 
ergy of £'sN — 10^^ erg, an ambient number density of 
n — 1 cm~^, and an acceleration efficiency of 6* — 0.1. 
The H.E.S.S. upper limits fall within the same order of 
magnitude of this value. 

SNRs with an upper limit on the p parameter above 0.1 
can be divided into two groups: those with low exposure 
and those with detection in VHE 7 rays. In the first case 
the H.E.S.S. coverage of the region is not deep enough 
to put a constraining upper limit on the SNR, whereas 
in the latter case VHE 7-ray emission has been detected, 
which can originate from the shell itself (e.g. G347.3-0.5, 
HESS J1713-397), from a possible bright pulsar wind neb- 
ulae (PWN, e.g. G0.9H-0.1, HESS J1747-281) or can re- 
sult from an overlapping bright source (e.g. G27.4H-0.0, 
HESS J1841-055). 

There are a few SNRs that appear to be underluminous, 
and they are interesting because, in principle, they have the 
potential of falsifying the model of hadronic 7-ray produc- 
tion. But the uncertainties in the input parameters to the 
model are large. In particular, the ambient density is usu- 
ally uncertain by an order of magnitude or more and the 
acceleration efficiency is also not well known. 

In summary, the present study constitutes a first falsifica- 
tion test for a commonly used model of the 7-ray visibility 
of supernova remnants. Further observations and future in- 
struments like the Cherenkov Telescope Array (CTA) with 
its vastly increased sensitivity and good survey capabili- 
ties will contribute to a better understanding of acceleration 
mechanisms of cosmic rays. 
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Table 1 : Parameters used for the SNRs considered within the present study, together with the preliminary values of the 
H.E.S.S. integral flux upper limit (at the 95% confidence level in units of the flux of the Crab nebula) and the derived 
upper limit on B ■ i?sN • n. The H.E.S.S. identifier is given for those sources where the SNR or an associated PWN is the 
most likely counterpart. 




